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We study scaling properties of the honeycomb fully packed loop ensemble associated with a 
lozenge tiling model of rough surface, when the latter is driven out of equilibrium by Kardar-Parisi- 
Zhang (KPZ) type dynamics. We show numerically that conformal invariance and signatures of 
critical percolation appear in the stationary KPZ state. In terms of the two-component Coulomb gas 
description of the Edwards-Wilkinson stationary state, our Hnding is understood as the invariance of 
one component under the effect of the non-linear KPZ term. On the other hand, we show a breaking 
of conformal invariance when the level lines of the other component are considered. 


I. INTRODUCTION 


Non-equilibrium stationary states in driven systems of¬ 
ten display scale invariance, and an important question 
is whether it can be extended to conformal invariance. 
In particular, on the 2-d plane, the existence of confor¬ 
mal invariance, which is infinite-dimensional, would have 
powerful implications on our theoretical understanding 
of non-equilibrium states of (2 -|- l)-d systems. Traces of 
conformal invariance are rare to find, and the answer to 
this question is unclear from a theoretical viewpoint. 

One remarkable positive result was obtained by Ber¬ 
nard et al [I| in the context of 2-d incompressible Navier- 
Stokes turbulence. Strong numerical evidence suggested 
that zero-vorticity lines behave as cluster frontiers of cri¬ 
tical percolation. The unexpected appearance of perco¬ 
lation in this context has raised great interest, yet still 
awaits a thorough understanding. In particular, does the 
presence of percolation teach us something on the turbu¬ 
lent state ? Is there a link between it and other contexts 
where critical percolation signatures occur, e.g^ nodal 
domains of wave functions in quantum chaos [2, 0 or 
kinetic Ising ferromagnet ? 

Recently, Saberi et. al. applied the idea (put forward 
initially by Kondev and Henley @) of studying loop en¬ 
sembles to the (2-|-I)-d Kardar-Parisi-Zhang(KPZ) equa¬ 
tion 0, 


dth = 


Vh 


+ 77 - 


( 1 ) 


This is a fundamental model of irreversible surface 
growth, and it also describes the turbulence without pres¬ 
sure IT], Isl , among other phenomena (for a recent review, 
see Q). Saberi et. al. observed numerically that level lines 
of stationary KPZ surface behave like self-avoiding walks, 
and so enjoy conformal invariance. This claim remains 
controversial and worrisome caveats [llj have not 
been resolved. 

Nevertheless, in this Letter we will show that, on KPZ 
rough surface models one can define a different loop en¬ 
semble which is conformal invariant and displays proper¬ 
ties of critical percolation. In particular, it passes strin¬ 
gent tests involving 3-point correlations, whereas KPZ 
level lines fail them. Moreover, we can understand our 
finding in terms of field theory. 



Figure 1: Left : A local snapshot of the unit cube/lozenge 
rough surface model, with the surface height h value on each 
vertex and level lines. Right : The fully-packed loops construc¬ 
ted on top of the configuration. 


II. MODEL AND NUMERICAL 
IMPLEMENTATION 


We will focus on the model whose configurations are 
tiling of unit cubes of the cubic lattice (or, in a 2d view¬ 
point, lozenge tiling of the plane). As shown in Fig. [1] 
left, a height function h is naturally defined on it. By lin¬ 
king middle points of each lozenge in this way (and 
similarly for the other two lozenge kinds, as in Fig. [1] 
right), one establishes a bijection between lozenge tiling 
configurations and honeycomb-lattice fully packed loop 
configurations. Olejarz et. al. proposed the following 
stochastic unit cube deposition process : for each time 
interval dt, an eligible move as shown below is chosen 

randomly : $ —>■ $. This dynamics is known to be in 
the KPZ universality class [Q : The coarse-grained li¬ 
mit of h satisfies Eq. [T] in the continuum. Symmetrising 
the above dynamics, we obtain a growth model descri¬ 
bed by the Edwards-Wilkinson(EW) equation,be., KPZ 

equation ([T]) with A = 0 : $ 

It is interesting to represent lozenge tiling configura¬ 
tions in terms of Ising spin configurations on the trian¬ 
gular lattice. In fact, see Eig. [D middle and right, gi¬ 
ven such a spin configuration which is maximally anti¬ 
ferromagnetic (be., for each elemental triangle, the spins 
on its 3 vertices are not the same), a lozenge tiling can be 
constructed by removing (from the triangular lattice) all 
the lattice edges connecting equal spins. The fully packed 
loops are then identical to the spin interfaces separating 
clusters of same spins. Therefore, standard methods for 
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constructing spin clusters and boundaries can be applied 
to both lozenge tiling models and triangular-lattice site 
percolation. (In terms of the height h, the Ising spin is 
obtained by a := exp(i7r/i), i.e., ±1 according to parity 
of h.) 

With this representation, the dynamics EW and KPZ 
are single-flip dynamics of Ising spins. In each trial, a site 
is picked at random and its spin cr is compared to that 
of its 6 nearest neighbours, (Tq, ..., us (in circular order, 
with (To pointing up). For EW, the central spin is flip¬ 
ped if (T = (To = (T2 = (74 or if (T = (Ti = (To = (75 (this 
corresponds to a Metropolis Monte Carlo dynamics); for 
KPZ, the central spin is flipped only if cr = (Ti = (T 3 = <75 
(this is a non-equilibrium dynamics). It is not hard to 
see that in terms of lozenge tiling, the dynamics just des¬ 
cribed is equivalent to the respective ones defined above. 
The number of lozenges of each type are preserved by 
the dynamics, and we choose to stay in the sector with 
equal number of lozenges of each type. This is assured by 
starting from the flat initial condition depicted in Fig. [3] 
middle. We let the system attain the stationary regime 
before constructing loop ensembles; snapshots of typical 
long loops are shown in Fig.[3]left(EW) and right(KPZ). 

Some care is needed to ensure the compatibility with 
the toroidal periodic boundary condition (p.b.c). Indeed, 
see Fig. [2] left, let ei and 62 be basis vectors of the trian¬ 
gular lattice edges forming 120 degrees(in terms of Eu¬ 
clidean inner products (ei,ei) = ( 62 , 62 ) = 1 ,( 61 , 62 ) = 
— 1/2), then the triangular lattice is partitioned into three 
sub-lattices generated by Ei = ei -|- 2 e 2 and E 2 = 61 — 62 , 
and the flat configuration is the one with up spin on one 
sub-lattice and down on the other two (see Fig.[3]middle). 
We fix our p.b.c for a system of size L by identifying 
LEi = LE 2 = 0, preserving the tri-partition and thus 
being compatible with the flat configuration. 

Data presented below are obtained on systems of size 
up to L = 2048. For each measure, le3 ^ led samples 
are generated, extending over a time scale corresponding 
to led ~ leS elemental moves per lattice site. 


A. Fully packed loops and level lines 

The stationary state of EW dynamics has the same 
weight for all tiling configurations and coincides there¬ 
fore with the well-known fully-packed loop model (which 
we shall call the EW FPL henceforth). Many results 
were obtained by Bethe Ansatz, [l^ and Coulomb gas 
approaches M- The EW FPL loop ensemble is confor¬ 
mal invariant and strictly related to critical percolation. 
In particular, the loops have the same fractal dimension 
Df = 7/4 (we define this notion below) as percolation 
frontiers. On the other hand, at the KPZ stationary state, 
the tiling configurations have different weights and no 
exact result is known. We will refer to the corresponding 
loop ensemble as KPZ FPL. We address the question 
whether the KPZ non-linear term affect the universality 
of EW FPL. For instance will the KPZ loop fractal di¬ 


mension be different from the EW one ? 

Let us recall that this is the case for level lines of /i, 
which are obtained by connecting lozenge middle points 
in the other way , see Fig. [T] left. In the EW case, 
these lines are described in the continuous limit by the 
level lines of Gaussian free field and have therefore fractal 
dimension D^''- =3/2 [T^. In the KPZ case, the fractal 
dimension of the same lines was estimated numerically to 
be Df''- w 1.33 i- Here we show that the critical perco¬ 
lation nature of the fully packed loop ensemble persists 
under the KPZ dynamics. 

III. SCALING EXPONENTS 

The loop fractal dimension is related to a scaling rela¬ 
tion between r, the distance between two extremities of 
a loop segment and I its length : 

(^ 2 ) _ 12/D,_ ( 2 ) 

This scaling behaviour holds when a I L, where 
a is the lattice spacing and L is the loop total length. 
In practice, we average over loop segments of length I < 
O.IL for L > 563. The result, shown in Fig. U main, 
suggests convincingly that fully-packed loops have D{ = 
1 jA: for both EW and KPZ stationary state. 

With the same loop segments at disposal, we also mea¬ 
sure winding angle variance. For a length-parametrised 
curve segment rs(s G [0,/]) of length Z, let 9g be a conti¬ 
nuous function satisfying dr/ds = (cos 0s, sin 0s), i.e., its 
tangent direction angle, then its (end-to-end) winding 
angle is defined as w := 0; — 0o. It is predicted that([l3|) 

{{0i-e,f) = j^^^\og{l/lo) (3) 

(where Zq is a lattice-dependent constant), holds for 
conformal invariant curves of fractal dimension Df. As 
shown in Fig. |T]inset, the above formula, with Df = 7/4, 
agrees perfectly with KPZ and EW FPL. This strongly 
supports conformal invariance in the KPZ FPL ensemble. 

We consider next the loop length distribution, expec¬ 
ted to have a power law behaviour 

(n(L))^L--+i. (4) 

The above scaling should be valid in an interval a <C L ^ 
Fmax, Fmax being the system size cut-off. The definition 
of the exponent r is the same as in @ . For critical perco¬ 
lation frontiers, it is known exactly that r = 15/7. In Fig. 
[S]main we superpose {n{L)) of KPZ and EW surface. The 
two distributions are indistinguishable in the scaling re¬ 
gime, where the power law is again in agreement with the 
percolation prediction. Using general scaling arguments, 
Kondev and Henley @ derived a relation between the 
roughness exponent of a scaling invariant surface, and 
the exponents Df and r associated to its level lines 

Dfir - I) = 2 - a. 


( 5 ) 
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Figure 2: Left : Generators of the trigular lattice ei,e 2 and that of a sublattice Ei,E 2 . Middle/right : A maximally anti¬ 
ferromagnetic spin configuration and its corresponding lozenge tiling. 



Figure 3: Left ■. Snapshot of some full packed loops in the stationary regime of EW. Middle : The fiat configuration. Right : 
Snapshot of some full packed loops in the stationary regime of KPZ. The snapshots are taken in system of size L = 64 ; distinct 
loops are drawn with different colors. 




Figure 4: Main : The mean square distance between extremi¬ 
ties (r^)(j/-axis) as function of loop segment length I (a;-axis). 
Inset : Semi-log plot of the winding angle variance {w^) as 
function of 1. Black lines are theoretical predictions 0, m 
with Df = 7/4:. 


Setting Df = 7/4, r = 15/7 one gets a = 0, which is 
true for any equilibrium critical loop model. But the fact 
that the same happens to KPZ is a priori unexpected. 
Note that this is not in contradiction with the non-zero 
roughness of the KPZ surface, as we recall that surface 
level lines form a different ensemble from FPL. One can 


Figure 5: Main : The distribution of loop length : n{L), the 
number of loops with length £ [L,L{1 -|- e))( 2 /-axis,e = 0.06) 
as a function of logL(a:-axis). Black Line is the theoretical 
prediction [4] with r = 15/7. Up inset : Roughness of the hid¬ 
den surface /i 2 : The mean squared difference of the hidden 
surface A/i 2 := ^ 2 (r) — /i 2 ( 0 ) in function of spatial separa¬ 
tion r = |r|. Down inset : Zoom-in of the microscopic regime 
where EW and KPZ FPLs are different. 


give a meaning to the exponent a in Eq. dSD by construc¬ 
ting a surface whose level lines are the FPLs. In practice, 
we assign a height increment crossing each loop by choo¬ 
sing randomly one of the values ±1. The obtained height 
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Table I: From left to right : a configuration contributing to 
Pi{xi,X2,X3), a configuration contributing to Pc(®i, ®2, xs), 
and a configuration not contributing to Pc{xi, X 2 , X 3 ) 

function [ 2 ^ we denote / 12 . This field is the second com¬ 
ponent of the vector-valued height function (/ii,/i 2 ) in 
the Coulomb gas [l^ approach to the FPL model. The 
other component hi is the lozenge height function h (see 
Fig. [T] left) introduced from the beginning. 

From a CFT viewpoint, Dotsenko et al [l^ remarked 
that EW FPL can be described by the following effective 
action : 

S[hi,h2]=Si[hi]+S2[h2], ( 6 ) 

where the hi and ^2 degrees of freedom are decoupled 
(the decoupling was first found numerically by El)- The 
action S 2 [h 2 ] describes the critical percolation (dense 
0{n = 1) loop model, 0), while is a free-boson 

CFT with central charge c = 1. The KPZ stationary 
state eludes any field-theoretical description, being far 
from equilibrium. Nevertheless, we put forward the hy¬ 
pothesis that correlation functions in the h 2 sector of the 
EW FPL are not affected by the KPZ dynamics, while 
those involving the field hi are. The results shown above 
supports this scenario. In particular, the field h 2 has lo¬ 
garithmic roughness for KPZ as well as for EW. This 
is also confirmed by explicit numerical measure of mean 
squared height difference in function of spatial separa¬ 
tion Fig. [S] up inset. In the following, we provide further 
evidence of the above hypothesis by looking at 3-point 
correlation functions. 


IV. 3-POINT STRUCTURE CONSTANTS AND 
CONFORMAL INVARIANCE 

Until now, our hypothesis is corroborated by obser¬ 
vables that can be related to 2-point functions. Yet truly 
profound predictions of CFT begin with 3-point func¬ 
tions. As illustrated in Table HI for /c > 1, we consi¬ 
der the probability Pi{xi,X2, ■ ■ ■ ,Xk) (/oop A:-point func¬ 
tion) that xi,X2, . ■., and Xk lie on the same loop, and 
Pc{xi,X2, ■ ■ ■ ,Xk) (cluster fc-point function) the probabi¬ 
lity that, for any couple of points Xi,Xj, there exists a 
curve connecting them and not crossing any loop. 

When k = 2, the CFT predicts Pn( 0 ,x) ^ \x\~‘^^° , o = 
c,l, with xi = 1/4 and Xc = 5/48 For k = 3 , We 



Figure 6: Main frames : 3-point ratios where the points x, y 
and 2 form a eqnilateral triangle of side r. Up : Loop 3-point 
ratio estimate Ci = 0.952(3). Measures using different mi¬ 
croscopic definitions of point are superposed. Down : Cluster 
3-point ratios compared to prediction Cc ~ 0.722 .... Insets : 
the corresponding 2-point function in function of distance r. 

consider the 3-point ratio 

r. _ Po{x,y,z) _, 

0(5 — , -? o — I. (') 

VPo{x, y)Po{x, y)Po{x, y) 

Global conformal invariance implies that Co is a constant, 
universal, which defines the CFT. 

For critical percolation, the exact value of Cc has been 
predicted by Delfino and Viti [2l| to be given by the Liou- 
ville CF T I22II . This conjecture was numerically confir¬ 
med [U, [ 2 ^. The exact value of C/ is not known (see, 
however, [2^), and to our knowledge numerical values 
are not published. 0 

We determine the ratios Cc and Ci for EW and KPZ 
FPLs. The 3-point ratios are calculated for x,y and 2 
forming equilateral triangles of varying radius; different 
geometries have been checked for the independence of Co- 
The results are shown in the main plots of Fig. [6l and 
compared to those obtained from critical site percolation 
frontiers on the triangular lattice. The numerics support 
strongly the following : 

• 2-point functions P(i,c) of EW and KPZ FPLs agree 
with the percolation exponents Xi = 1/4, ccc = 5/48. 

• It exists a scaling regime where the 3-point ratios are 
constant for EW and KPZ FPLs. 
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• The value of Cc, which is the same for EW and KPZ 
FPLs, and the numerical result is in excellent agree¬ 
ment with Delfino-Viti’s prediction for the critical 
percolation : 

Cc = 0.722 ... (8) 

• The loop structure constant Ci is the same for critical 
percolation, EW and KPZ FPLs, its numerical value 
is estimated to 

Cz= 0.952(3). (9) 

We further collapsed this estimate against an inde¬ 
pendent measure on critical bond percolation on the 
square lattice. 

The above results are in agreement with our hypothe¬ 
sis that the conformal invariance in the /12 sector is not 
broken by KPZ dynamics. Nonetheless, we remark that 
the EW-KPZ coincidence of the observables considered 
occurs not only in universal exponents, but often also in 
non-universal amplitudes. This appears far from trivial 
and begs an elementary explanation. 

V. hi SECTOR OBSERVABLES 

Contrary to what can be suggested by the results of the 
previous section, we emphasize that the EW and KPZ 
FPLs are not identical. It indeed suffices to encode some 
information on hi in the FPLs to see distinct behaviours 
for EW and KPZ. Here we show how to construct such 
observables of individual loops. For this, consider a loop 
segment of length /. It can be described by a sequence of 
left/right turns of 60 degrees, d0s := 0^ — 0^-1 = if j s = 
1,... ,1. (The fully packed loops live in fact on the hexa¬ 
gonal lattice dual to the triangular lattice.) Recall that 
the winding angle, defined as associated to the 

loop segments cannot distinguish KPZ from EW for I 
large. However, the observable 6h{l) := f 
does so, see Fig. [7] left. The variance {{Sh)'^) has loga¬ 
rithmic growth for EW and power-law for KPZ. Indeed, 
it can be shown that |5/i(^)| = |/ii(ro) — hi{ri)\ where 
vq and ri are the extremity positions of the loop seg¬ 
ment. Thus, 6h{l) encodes the roughness of hi (logarith¬ 
mic for EW, positive power law for KPZ) into fully pa¬ 
cked loops, explaining the observed distinct scaling be¬ 
haviour. Note that the definition of Sh differs from the 
winding angle w just by an alternative sign. To produce a 
second example, we consider the extremity displacement 
^(0 •= Sl=i where dr^ = (cos0s, sin0s), which can¬ 
not distinguish KPZ and EW at long distance. On the 
contrary, the quantity Sd{l) := grows with 

distinct power laws for KPZ and EW fully packed loops 


(Fig. [7] right). Again, one can realise that this measure 
encodes information from hi, whose roughness is also 
reflected in correlation of lozenge orientation. Denoting 
no, ni and 712 the number of lozenges of each orientation 
passed by the loop segment, it is not hard to see that 
(Sd)"^ = |no-I--I-n 2 e“*^/^’^| and measures the 
deviation of nz’s from the expected average 1/3. 

We finally consider the 3-point ratio for the level lines 
(of height hi), which is of more physical interest, see Fig. 
[T]left. Saberi et. al. conjectured that the KPZ level lines 
are in the universality class of self-avoiding walks. Here 
we compare loop 3-point ratio of KPZ level lines with 
that of external perimeters of critical site percolation 
clusters, since they are believed to be identical to self¬ 
avoiding walks As we show in Fig. [51 percolation 

external perimeters display a well-defined loop structure 
constant, whereas KPZ level lines do not. Finally, we ob¬ 
serve that the loop 3-point ratio for EW level lines attains 
also a constant. This result, taken together with another 
incompatibility (noticed in 0 ) with hyper-scaling re¬ 
lations in and estimates of KPZ roughness exponents 
[13, casts doubt on the conformal invariance of KPZ level 
lines and their conjectured identification to self-avoiding 
walks. 


VI. CONCLUSION 

In this work, we have studied KPZ stationary state 
compared to the EW one on a natural discrete model. 
We have shown that, in terms of the two-component sur¬ 
face description, whereas hi undergoes dramatic change 
under the KPZ dynamics, the “hidden” component h 2 re¬ 
mains invariant. The fully packed loop ensemble associa¬ 
ted with h 2 enjoys conformal invariance, and has critical 
properties indistinguishable from that of critical perco¬ 
lation. On the other hand, conformal invariance appears 
broken for level lines of hi. The decoupling of /12 from hi 
in KPZ occurs exactly in the scaling regime, even if at the 
microscopic level, there is no a priori reason. We are left 
to wonder what implications, if any, our observations can 
have on field-theory of non-equilibrium stationary states 
in general, and on the genuine KPZ phenomenology in 
particular. 
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Figure 7: Left : the observable {3h)‘^ as function of loop segment length 1. Right : the observable {3d)^ as function of loop 
segment length 1 . 



Figure 8: Loop 3-point ratios of EW and KPZ level line en¬ 
semble as well as the external perimeters ensemble of critical 
percolation. Different curves of the same type indicate the 
error. 
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